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Abstract. One more mode developed to get eigen energies and states for the 
one-electron Dirac's equation with spherically symmetric bound potential. For the 
particular case of the Coulomb potential it was shown that the method is free of so 
called spurious states. The procedure could be adapted to receive highly exited states 
with great precision. 



1. Introduction 

It's difficult to overestimate the importance for the relativistic calculations numerical 
approaches to get full spectra of Dirac's equation. One of the most successful is so 
called B-spline approach (Carl de Boor 1978), (Johnson et al. 1988). While numerical 
implementation of the problem so called spurious states can arise. Those states need 
special treatment (Shabaev et al. 2004). All said above prompts to explore new methods 
to solve Dirac's equation numerically. Namely in this paper we'll investigate the 
advantages that squared Hamiltonian is able to give for reaching the goal. Throughout 
in the paper atomic units system {h = e = m) is utilized. 

2. Squared Dirac's Hamiltonian 

Let's take the squared one-electron Dirac's Hamiltonian in the form: 



Where h = c{a ■ p) + V + mc^ (3 is a well-known representation of the one-electron Dirac's 
Hamiltonian (Labzowsky et al. 1993) with some bound potential V and e is an arbitrary 
real-valued parameter. Having deployed the ([1]) everybody can easy receive 




(1) 



H = (?{a ■ -pf m^c^ + + mc'{{a ■ p), /?} -h c{(a ■ p), V}^ 
2mc^pV - 2eh + 



(2) 
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(3) 



where {a, b} = ab + ba is a positive commutator of the operators a and b. Keeping 
in mind that the fourth term in the above formulae ([2]) amounts to zero (due to the 
properties of the Dirac's matrixes) one can finally get 

H = -ffc^ A +m2c^ + V^ + 2mc^/3V - ihc{{a ■ V),V}+ 
2iche{a ■ V) - 2eV - 2emc^(3 + e 

where for the Coulomb potential V = — e and Z are an electron charge and charge 
number of the nucleus respectively. Our further purpose is to solve the eigen problem 
for the squared Dirac's equation 

= (4) 

numerically. It's clear that X = [E — e)^ where E is an eigen number for the h operator. 
In order to reach that goal one can adapt the usual representation for the eigen function 
of the Dirac's equation with spherically symmetric bound potential (Labzowsky et 
al. 1993). Namely let's take 

that is well-known representation for the Dirac's bispinor, g = g{r) and / = /(r) are 
upper and lower radial component functions respectively. For the purpose of numerical 
implementation one needs to expand those functions over finite set of basis functions. 
We used here finite basis constructed from the B-splines (see for example (Johnson et 
al. 1988)). So if one has got the finite basis set {Bi{x)}^^i it's straight forward to gain 
the generalized symmetric eigen value matrix problem for the equation (j4]) 

Where H^^\ H^^\ B^-^\ are the sub-matrixes {n x n) of the general 

matrixes H and B {2n x 2n), with the matrix elements 

4^ = 4^ = (B^r^B,) (7) 

= -h'c' {{B.r'B'D + 2{B,rB'^) - 1(1 + 1){B,B^)) + 

m^c'^iBir^Bj) + (Bir^V^Bj) + 2mc'^{Bir^VBjy (8) 

2emc^{Bir^Bj) - 2e{B,yVB^) + e'^{Bir^Bj) 

H\f = -ch {2{Bir^VB'.) + 2(1 - K){BirVB,) + {Bir^V'B^)) + 
2ech{{Bir^B'^) + {l-K){B,rBj)) 

H\f = ch {2{ByVB'.) + 2(1 + K){B,rVBj) + {Bir^V'B,)) - 
2ech {{Bir^B'j) + (1 + K)iBirB,)) 



(9) 
(10) 



Hjf = -h'c^ {{B^B';) + 2{B,rB'^-l(l+l){B,B^)) + 

m^c'^iBir^Bj) + (ByV^Bj) - 2mc^{ByVB-)+ (11) 
2emc'^{Bir^Bj) - 2e{Bir^VBj) + e\Bir^Bj) 



besides (/) = /„°° f{r)dr. 
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3. Spectra quality test, e = operator 

To examine the grade of the spectra Bi H-like ion was chosen. The knot net sequence for 
the B-spIine construction was taken in accordance with the formula k(i) = i^ + a. 
The boundary conditions were adapted and /(a) = f{h) = g{a) = g{h) = in order 
to prevent the kinetic energy not to be negative. The outcome of the calculation 
(diagonalization of (jl]) with e = ) is presented in Table [H In Table [1] the lowest 
eigen values are compared not only with the corresponding eigen values gained from 
Sommerfeld formula but also with the eigen values received from the virial relations 
( |T2| [T3l) . The virial theorem for the operator h and yields formulae 

E = mc^{^\(3\-<^) (12) 

and 

E"^ = m^c^ + mc^{m\V(3\^!) (13) 

respectively. As one can see from Table [H the agreement is quite good. In addition to 
this different spectra (especially with k > 0) were analyzed and it was found out the 
absence of so called spurious states. To avoid problems with spurious states special care 
needs while numerical implementation of the equation = E'^ (Shabaev et al. 2004). 



Table 1. Several lowest energy levels for Bi H-like ion with I = 0, j = ^. The order 
of splines = 6 and the number of intervals for the B-splines construction is n = 130. 
Parameters for the knot sequence construction are a = 10^^^ and 6 = 10 respectively. 
E stands for self energy value gained by diagonalization of the equation ^ . Eg means 
values received from Sommerfeld's formula. The virial theorem tests are also presented. 
The Eyri and Eyr2 columns give the self energies received with formulae ([T^ and (fTS]) 
respectively. 



Es 


E 


E 1 


Evr2 


-3836.36956 


-3836.36901 


-3836.36498 


-3836.36583 


-984.921152 


-984.921036 


-984.920525 


-984.920581 


-425.768940 


-425.768909 


-425.768766 


-425.768777 


-234.607639 


-234.607627 


-234.607570 


-234.607573 


-148.000373 


-148.000368 


-148.000340 


-148.000340 


-101.712322 


-101.712319 


-101.712304 


-101.712304 


-74.1444640 


-74.1444625 


-74.1444530 


-74.1444528 


-56.4227560 


-56.4227552 


-56.4227489 


-56.4227487 


-44.3652984 


-44.3652978 


-44.3652936 


-44.3652933 


-35.7940987 


-35.7940983 


-35.7940953 


-35.7940951 


-29.4849330 


-29.4849328 


-29.4849305 


-29.4849304 
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4. The states from continuous part of spectra 

As one can see the outcome of the diagonahzation of the equation (jlj) with e = 
yields for the eigen values corresponding to the continuous part of spectra to be doubly 
degenerated. In general case the wave function \E' (for that part of spectra) is a mix 
of the couple of functions y?+ (really the self energy function of the h operator with 
the corresponding self energy E > 0) and (/?_ (really the self energy function of the h 
operator with the corresponding self energy —E) 

^ = Ciip+ + C2ip^ (14) 

with some arbitrary coefficients Ci and C2. Utilizing the properties of the self energy 
functions (fTSl). (fT6l) and (ITTl) 

h(p+ = E(p+ (15) 

/iV5- = -E^_ (16) 

/i^ = CiEip+ - C2Eip^ (17) 

one can easy separate the positive and negative spectra functions. 

/i^ + S$ = 2CiE(p+ (18) 

- E"^ = -2C2^V5- (19) 

Finally the formulae ( |T8l) and ( |T9l) present the desired functions. Everyone can also easy 
see that those functions are orthogonal (They really must be!). The formula (!20|) gives 
the proof of the fact. 

{ip^\ip_) = -{ACiC2E^)-^{h'^ + - = 

-{4CiC2E^)-\'^\{h + E){h - S)|^) = (20) 
-{ACiC2E^y^{'^\h^ -E^\^) = 

To conclude this section we can say that there is no problem to get full spectra and 
eigen functions of Dirac's one-electron Hamiltonian from the squared one. 



5. The e 7^ operator 

The case of the operator (jlj) with e 7^ is very interesting one. When e = the 
spectra of the operator (j4j) looks like one of the Schrodinger equation, namely we have 
two part spectra, the lowest energies are bound states and energies of the states from 
continuous spectra above them. The e 7^ transforms spectra of Dirac's equation even 
more, namely appropriate choice gives possibility to make any energy level the lowest 
one. This possibility was tested and gave positive results. The most important thing is 
that eigen functions and values (energies) could be found separately from each other as 
the lowest eigen values of the squared Hamiltonian (jlj). 



Another method to solve Dirac's one-electron equation numerically 



5 



6. Conclusions and future perspective 

The results of the present work are following. Another way to get spectra and 
corresponding eigen functions of one electron Dirac's equation is presented. It's shown 
the absence of spurious states. It was pointed the very attractive method to get highly 
excited states of the Hamiltonian separatly each other as the lowest eigen values of some 
squared Hamiltonians. The further purpose of author is a generalization of the method 
for many electron case. 
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